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Abstract. We consider a system of parallel straight edge dislocations and we analyse its as- 
ymptotic behaviour in the limit of many dislocations. The dislocations are represented by points 
in a plane, and they are arranged in vertical walls; each wall is free to move in the horizontal 
direction. The system is described by a discrete energy depending on the one-dimensional hori- 
zontal positions x; > of the n walls; the energy contains contributions from repulsive pairwise 
interactions between all walls, a global shear stress forcing the walls to the left, and a pinned 
wall at x = that prevents the walls from leaving through the left boundary. 

We study the behaviour of the energy as the number n of walls tends to infinity, and charac- 
terise this behaviour in terms of T-convergence. There are five different cases, depending on the 
asymptotic behaviour of the single dimcnsionlcss parameter f3„, corresponding to /3„ < 1/n, 
1/n <C fin <C 1, and n S> 1, and the two critical regimes /3 n ~ 1/n and n ~ 1. As a con- 
sequence we obtain characterisations of the limiting behaviour of stationary states in each of 
these five regimes. 

The results shed new light on the open problem of upscaling large numbers of dislocations. 
We show how various existing upscaled models arise as special cases of the theorems of this 
paper. The wide variety of behaviour suggests that upscaled models should incorporate more 
information than just dislocation densities. This additional information is encoded in the limit 
of the dimcnsionlcss parameter /? n . 



1. Introduction 

1.1. Dislocation plasticity. One of the hard open problems in engineering is the upscaling of 
large numbers of dislocations. Dislocations are defects in the crystal lattice of a metal, and their 
collective motion gives rise to macroscopic permanent or plastic deformation. 

For systems of millimeter-size or larger there is a fairly complete theory of macroscopic plasticity, 
in which dislocations are not modelled cxplictly (see e.g. [27, 29, 28, 5]). For smaller systems, 
however, the so-called size effect [25, 34, 20] suggests that it is necessary to take the distribution 
of dislocations into account. In this point of view the size effect arises when the length scale of 
the system becomes similar to the typical scale at which the dislocation density varies. 

To address these small-scale effects a number of competing (mainly phenomenological) dislo- 
cation density models have been derived by upscaling large numbers of dislocations (e.g. [16, 17, 
21, 22, 23, 31, 35]). The unknowns in this type of model are various types of dislocation densities, 
whose evolution in time is described via conservation laws equipped with constitutive laws both 
for the velocity of the dislocations and for their interaction. 

The use of densities (as opposed to keeping track of the behaviour of each dislocation) seems 
reasonable, since the typical number of dislocations in a portion of metal is huge. For topological 
reasons dislocations are curves in three-dimensional space, and therefore the density of dislocations 
has dimensions of length/volume or m~ 2 . A dislocation density of 10 15 m~ 2 (typical for cold-rolled 
metal [30, p. 20]) translates into 1000 km of dislocation curve in a cubic millimeter of metal. This 
high number explains the interest in avoiding the description of the individual behaviour of each 
dislocation, and focussing on the collective behaviour instead. It also explains the general belief 
that this should be possible. 
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The research done here, however, suggests that the situation is more subtle. It was triggered by 
the earlier study [36] . The outcome of [36] and the results in this paper suggest that the dislocation 
density alone is not capable of describing the evolution of large numbers of dislocations. To put 
it succinctly, a density simply does not contain enough information to characterise the behaviour 
of the system, even in aggregate form. This is because the density only characterizes the local 
number of dislocations per unit area, and needs to be supplemented with more information on 
their spatial arrangement in order to give a satisfactory answer. We show below how this point 
arises from the results of this paper. 

A separate reason for the analysis of this paper is the uncommon form of the energy. Although 
it is a simple two-point interaction energy in one spatial dimension, the behaviour in the many- 
dislocation limit does not fit into any of the standard cases as described e.g. in [12]. This is due 
to the combination of all-neighbours-interaction (each pair interaction is counted, regardless of 
distance), and an interaction potential that is globally repulsive. 

1.2. The model of this paper. We consider a system of pure edge dislocations whose dislocation 
lines are straight and parallel to one another as in [36]. These dislocations can be modelled as 
points in the plane orthogonal to the direction of the dislocation lines, and this identification 
has been done systematically in the literature. The slip planes are horizontal, i.e. parallel to the 
i-coordinate, which implies that the dislocations can only move in the horizontal direction (see 
Figure 1). In addition, the dislocations are organized in vertical walls with a uniform spacing 
of size h (in m). In the model below there will be a finite number n of such walls, which each 
will extend indefinitely in the vertical direction. The total degrees of freedom of the system are 
therefore the horizontal positions < X\ < • • • < x n (in m) of the walls. A constant global 
shear stress forces the walls towards a fixed barrier which is modelled as an infinite wall of pinned 
dislocations at xq = 0. 




Figure 1. The dislocation configuration considered in this paper. Infinite, ver- 
tical walls of equispaced dislocations are free to move in the horizontal direction. 
A wall of fixed dislocations is pinned at xq — and acts as a repellent. 

We assume that the dislocations are spaced significantly farther apart than the atomic lat- 
tice spacing, which implies that the interactions between dislocation walls are well described by 
conventional formulae based on linear elasticity. 

Given these assumptions, the system is driven by the discrete energy 

K n n f * ■ - v ■ \ n 

£(5i,.-.,in) = f EEM^V^ +ff L>' a- 1 ) 

8=1 j = \ ' »=1 

where K := Gbn/2(1 — v) 1 G [Pa] is the shear modulus, b [m] the length of the Burgers vector, 
v [1] the Poisson ratio of the material, and a [Pa] the imposed shear stress. The (dimensionless) 



DISLOCATION DENSITY 3 

interaction energy V is 

V{s) := I SC oth^- ^log(2sinh7r S ) = - J* 1 - -I log(l - e^), (1.2) 
and its derivative is the (dimensionless) force exerted by a wall on another wall at distance s, 
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FIGURE 2. The interaction energy V. 

The first term in the discrete energy £ in (1.1) is fully repulsive: each pair of walls repels 
each other, with a potential that diverges logarithmically as the walls approach each other (see 
Figure 2). The second term of the energy, accounting for the global shear stress, drives the walls 
to the left. The repelling nature of the left barrier is implemented by pinning a wall at Xq = 0. 

Stationary points of this energy are equilibria of the mechanical system, and under the assump- 
tion of a linear drag relation (see e.g. [30, Sec. 3.5]) the evolution of the system is a gradient flow 
of this energy. 

Although the model is highly idealised, it has a number of properties that make it both in- 
teresting and not unrealistic. The fact that multiple dislocations move along exactly the same 
slip plane is natural, because of the way they are generated from Frank-Read sources (e.g. [30, 
Sec. 8.6]). Moreover, although the assumption of an arrangement in equispaced vertical walls is 
clearly an idealisation, it is on the other hand not unrealistic since equispaced vertical walls are 
minimal-energy configurations. Walls of edge dislocations are locally stable, in the sense that if one 
of the dislocations deviates from its wall position, either horizontally or vertically, it experiences a 
restoring force from the other dislocations that pushes it back. Finally, the vertical organization 
in walls is also justified by correlation functions calculated from numerical simulations (e.g. [40]). 

Another interesting aspect of this model is that existing, phenomenological dislocation-density 
models can be applied to it to give predictions of the upscaled behaviour — which can then be 
tested against the rigorous results of this paper. In Section 1.5 we discuss three of these, whose 
predictions for this system are summarized in Table 1.2: 
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Reference 


Stationary state 


<7int 


Head & Louat 1955 [26] 


/>(*) - 4^ 




Groma, Csikor & Zaiser 2003 [23] 


p(x) = Ce~™ 


-d x p/p 


Evers, Brekelmans and Geers 2004 [19] 


p(x) — C — ax 


-d x p 



Table 1 . Various predictions for the limiting behaviour of this system of dislo- 
cations. The system is characterized by the density p of dislocations; a- mt is the 
prediction of the stress field generated by p. Parameters have been absorbed into 
C and a for simplicity. We give a full discussion in Section 1.5. 



As we show below, the results of this paper allow us to make sense of the three different 
predictions in this table. 

1.3. Main result. The main mathematical result of this paper is the characterization of the 
limit behaviour of £ as n — > oo. Since this behaviour depends strongly on the assumptions on the 
behaviour of the whole set of other parameters in this system, that is h, K (or G, 6, and u), and a, 
we assume that all parameters depend on n. In fact the parameter space is only one-dimensional, 
since the problem can be rescaled to depend only on the single dimcnsionlcss parameter 



Pn:=J^- = J 9 ,; Gnb ? . ■ (1.3) 
V na n h n y 2n(l - v n )o n h n 

In mechanical terms, /?„ measures the elastic properties of the medium (described by K n ) in 
comparison with the strength of the pile-up driving force u n . Large /?„, therefore, corresponds to 
weak forcing, and small [3 n to strong forcing 

We characterize the limiting behaviour of the system by proving five T-convergence results, for 
five regimes of behaviour of j3 n as n — > oo, after an appropriate rescaling of £ and (xi, . . . ,x n ) 
(rescalings that lead to the functional E n k \x\, . . . , x n ), for k <E {1, 2, 3, 4, 5}, defined in Theorem 
1.1). A consequence of T-convergence is the convergence of minimizers. 

Both the T-convergence of the energy and the convergence of minimizers depend on a concept 
of convergence for the set of wall positions (xi)f =1 , or their rescaled versions (xi)™ =1 , as n — > oo. 
A natural concept of convergence for such a system of wall positions is weak convergence of the 
corresponding empirical measures (which we prove being equivalent to the weak convergence of 
the linear interpolations of the wall positions in the space of functions with bounded variation, see 
Theorem 2.2). For a vector x € K™ define the empirical measure as 

1 ™ 

»n = -y^Sxi- (1-4) 

i=i 

The weak convergence of p n to p, written as p n — ^ p, is defined by 

/ if(y) p n (dy) "—^2? / tp(y) p(dy) for all continuous and bounded ip, 
Jn Jn 

where := [0, oo). This is the concept of convergence that we use in this paper. 

Theorem 1.1 (Asymptotic behaviour of £). In each of the cases below, boundedness of the func- 
tional En implies that the empirical measures p n :— ^ Y^i=i are compact in the weak topology. 
In addition, the functional E n k ^ T-converges to a functional E^ with respect to the same weak 
topology. 

Case 1: If (3 n <C 1/n, then define the rescaled positions X\,...,x n in terms of the physical 
positions xi,...,x n by 

x i — X{ (1-5) 
nK n 
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and define 

E^\xi,...x n ) := -5— £(x!, . . . ,x n ) + —^(log2irn 2 /3l - I). 



Then T-converges to 



e {1) {p) ■= J J lo § \ x ~v\ v(dy)K dx ) + J x v{dx). 

Cases (2—4): If (3 n <~ 1/n, 1/n <C /3 n <^ 1> 0J " Ai ~ 1, then define the reseated positions 
X\ , . . . , x n as 



^V-w (L6) 

and define 

^(pCli • • • j • / o , £(^1) ■ • • i x n)- 

yJn A K n h n a n 

Then En ^ T-converges to 
e{2) (p) - =< ^ v ( c ( x - V)) H{dx)n(dy) + / xfj,(dx) if n/3 n -> c, 

, ~ } ,«.r 

+ oo otherwise J 

(1.8) 



c / Vr.fi ( ) p{x)dx + / xp(x)dx if p,{dx) = p(x) dx 



E^(p):={ J n cn \p(x)J yy J J n PW JAn J, i//3„^c, 

+ oo otherwise 



(1.9) 



where the function V e s in (1.9) is defined as 

oo 

F cff (s) :=^F(fcs)- 



fc=i 



Case 5: If f3 n ^> 1, then define the reseated positions x\,. . . ,x n as 

Xi = Xi ( — nh n log ^ 11 

and define 



2tt V 7r nh n a n 



E^(x u ...x n ) := (^n 2 /^ log (f^)) ■ ■ ■ >*») 



(5) 

TTien £J„ T-converges to 



i i xp(x)dx if u(dx) — p(x) dx and p<l C-a.e. 

E {5 \p):=lJ n >, (1.10) 

[ + oo otherwise J 

where £ is the Lebesgue measure. 

The limiting energies have the nice property of strict convexity, either with respect to the linear 
structure in the space of measures, or in the sense of displacement convexity [32]. This gives 
uniqueness of minimizers: 

Theorem 1.2 (Existence and uniqueness of limiting minimizers). For each k e {1,2,3,4, 5} ; E^ 
has a unique minimizer in the set of non-negative, unit-mass Borel measures on fi. 
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As a consequence we can characterize the behaviour of sequences of minimizers: 

Corollary 1.3 (Convergence of discrete minimizers). Let the asymptotic behaviour of f3 n be as 
in case k £ {1,2,3,4,5} of Theorem 1.1. Let (x",...,x™) be a sequence of n-vectors such that 
for each n, (x",..., x™) is minimal for £. Then, reseating (x™,...,x™) to (x™, . . . , x") as in 
Theorem 1.1, the corresponding empirical measure \i n — j % J27=i ^™ converges weakly to the global 
minimizer of E^ . 

The proof of Theorem 1.1 is the subject of Section 3; Theorem 1.2 and Corollary 1.3 are proved 
in Section 4. 

Figures 3-7 show some numerical examples of the matching between discrete and continuous 
energies. Note that the optimal discrete density p n plotted in the Figures below is defined for 
every i = 2, . . . , n — 1 as 

i \ 2A n 

Pn(Xi) := , 

X^-j-l X^_i 

where (xi,...,x n ) is the minimiser of the discrete energy E„ fc , for k = 1,...,5 and A n is a 
normalization factor ensuring that the area below the linear interpolant of p n is one. 



500r 

450 - 

400 - 
350 - 
300 - 
250 - 

200,; 
150, r 

ioo| 

501 



0.05 0.1 0.15 0.2 0.25 

Dislocation positions 

FIGURE 3. Optimal densities relative to E^ and E^\ for n = 150 and {3 n = 6/(n^/n). 



1.4. Five regimes. The role of j3 n and of the different asymptotic regimes can be understood as 
follows. Define the average dimensional distance between two walls (assuming n even) as 

Ax:=% 
n/2 

Note that x n /2 is a 'middle' wall, and therefore a reasonable indication of the size of the pileup. 
Assuming cases 2-4, we can then rewrite (1.6) as 

Ai 

— — 2x n l2Pn- 

K 

If the empirical measures /i n in (1.4) converge, then x„/2 = 0(1); this equality therefore indicates 
that p n is a measure of the aspect ratio Ax//i„, or, put differently, nf3 n is a measure of the total 
length of the pileup, relative to h n . 

Cases 2-4, therefore, can be understood heuristically as follows: 
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Figure 5. Optimal densities relative to and for 
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If (3 n — > and nf3 n —> oo (case 3, and Figure 8(b) below), then the range of the ratio 
\xi — Xj\/h n , which appears as an argument of V in (1.1), asymptotically covers the whole 
range from to oo. In this case the discrete system effectively samples the integral J V, 
and this integral therefore appears in the limit energy (1.8). 

If /3„ — > c > (case 4 and Figure 8(c)), then the sampling of V does not refine, but 
remains discrete, and instead of the integral J V we find the discrete sampling V e g (1.9). 
If nf3 n — > c > (case 2 and Figure 8(a)), then the pile-up is not long enough to cover the 
whole of the integral of V. In addition, in this case the length scales of \x n and of V are 
exactly the same, and a convolution integral results. 
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FIGURE 6. Optimal densities relative to E n ' and E^ A \ with n = 150 and j3 n = 1 
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Case 1, where /3 n is so small that nj3 n — > 0, is a variant of case 2, but now the dislocation 
walls are pushed completely into the logarithmic singularity of V at the origin (since by (1.5) the 
typical total length of the pile-up is nK n /a n = h n n 2 [3 n , and this is small with respect to h n ). We 
observe that by the definition of f3 n (1.3) this situation corresponds to strong forcing, which pushes 
the dislocation walls closer to each other. Because of the scaling dependence of the logarithm, a 
multiplicative rescaling in space (in order to make the sequence \i n converge to a non-trivial limit) 
results in an additive rescaling of £ . The corresponding picture is similar to Figure 8(a). 

In case 5, where j3 n is large, the value of (xi — Xj)/h n also becomes large; even the two closest 
dislocation walls have distance asymptotically larger than h n . Then the dislocations only sample 
the exponential tail of V. By the definition of (3 n (1.3) this case corresponds to very weak forcing. 
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FIGURE 8. Cases 2-4 of Theorem 1.1. Since f3 n is a measure of the aspect ratio 
Ax/h n , the scaling of f3 n determines which values the ratio (£i — Xj)/h n takes in 
the argument of V in (1.1). 

The degenerate nature of the limit of the interaction energy, which is zero if /i = pdx with p < 1, 
and +oo otherwise, arises from the 'winner takes all' behaviour of the exponential function. 

Other possibilities for j3 n . One might wonder whether other scaling behaviour of /?„ could give 
different results. Although it is certainly possible to construct sequences /3 n that do not fit into 
the five classes above, by taking subsequences one can reduce the behaviour to one of these five 
possibilities. Of course, if different subsequences have different asymptotic behaviour, then one 
does not expect the functionals to converge; this non-convergence is a further indication that one 
should separate the cases by dividing into subsequences. 

1.5. Comparison with mesoscopic models in the engineering literature. As mentioned 
above, one motivation for this research is the derivation of a model describing the behaviour of 
densities of dislocations from a more fundamental microscopic model described by the discrete 
energy (1.2). The need for a rigorous derivation of such a dislocation-density model is underlined 
by the fact that multiple models exist in the literature (see Table 1.2) that are inconsistent with 
each other and whose range of validity is not clear. 

In the case of this paper, straight parallel edge dislocations in a single slip system, the upscaled 
evolution equation for the dislocation density (or measure) [i is expected to be of the form 

d t p + d x (vfi) = 0. (1.11) 

Here v(x,t) is the velocity of dislocations at (cc, t) , and is usually taken to be 

v=^-{<r- mt -a). (1.12) 

Here B is a mobility coefficient, a is the externally imposed shear stress (as above) and o~i n t is the 
shear stress field that the dislocations themselves generate, which is assumed to depend on the 
dislocation density and on the gradient of the density, i.e., a- lrA = cri nt (x, p, d x fj,). The structure 
(1.11)-(1.12) arises naturally from the evolution equations for the discrete system, 

= --d Zi e(*u • • ■ ,*») = - m £ - ^ (1.13) 

which suggests that a- ln t should be the upscaled limit of the interaction forces, represented by the 
sum in (1.13). 

The different models proposed in the engineering literature differ in the form of the internal 
stress (Ti n t they suggest, as shown by Table 1.2, and the arguments leading to the specific choice 
of (Tint can not always be rigorously justified. 
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In contrast to these mostly phenomenologically derived expressions, the convergence results of 
Theorem 1.1 offer a rigorous characterization of the limiting internal stress <7i nt in the different 
cases corresponding to our limit models k = 1, . . . , 5. The Euler-Lagrange equations of the limit 
functional E^ k \ expressed in terms of the measure (i or the density p, are (taking c = 1 for 
simplicity) : 

(k = 1) -\d x (log * M ) + 1 = 0; *£> = \d x (log * M ) ; 

(k = 2) d x (V * p) + 1 = 0; a£l = -d x (V*n); 

(k = 3) ^ jf V(t) dtj d xP +l = 0; = - Qf dt) x p; 



- 4 ) ^ ^ ( J) - 1 = 0; = 4j V" s Q ) ft ,, 



We leave out the case fc = 5 since its Euler-Lagrange equation is too degenerate to be useful. 

The Euler-Lagrange equation in case k = 1 coincides with the one derived by Eshelby, Frank, 
and Nabarro [18] and Head and Louat [26] in the case of n dislocations in one slip plane-rather 
than n dislocation walls. This is consistent with the fact that when fi n -C 1/n, the dislocation 
walls are much closer to each other horizontally than the vertical spacing h n (see the discussion 
of fi n above), and therefore an approximation by a single-slip-planc setup seems appropriate. 

The internal stress a^l coincides with the one proposed by Evers, Brekelmans and Geers [19]. 
As far as we know, the limiting energies for k = 2, 4 and the internal stress associated with 
them have not been mentioned in the engineering literature yet. 

Groma, Csikor, and Zaiser [23] derived the internal stress cr^ CZ ^ = —d x p/p (up to constants) 
starting from a discrete distribution of dislocations where the horizontal and vertical separation 
of the dislocations is of the same order. In our formulation this corresponds to the case k = 4, 
fi n ~ 1. Therefore it is interesting to compare a^^^ with o^. As it turns out, cr^ t CZ ^ can 
be formally obtained from a[^ t by making two approximations. The first consists in disregarding 
the interaction between walls that are not nearest neighbours, which is equivalent to replacing the 
effective potential V^ s with V'. The second approximation is to substitute the force V'(s) with its 
first-order Taylor-Laurent expansion close to zero, namely — Via these two approximations 

a^l reduces (up to a constant) to o-^ t CZ \ 

Although this derivation can formally be made, it can not be made rigorous, since the two 
approximations are mutually incompatible. In fact, if Taylor-expanding V is sensible, then the 
logarithmic singularity of V implies that there is interaction with all neighbours, as is always the 
case for logarithmic interactions (and as is the case for case k = 1 above). Therefore neglecting 
all but nearest neighbours is unjustified. Moreover, the Taylor expansion of V^ s and of V close 
to zero are quite different since when s is small, 1 V(s) w — 1/n 2 log \s\, while 



As] Jw 



Ves(s) « — / V. (1.14) 



Therefore truncating to nearest neigbours (replacing V e g by V) and then Taylor-expanding V 
results in a large error. We refer to the companion paper [37] for further discussions on this point 



^This follows from the two inequalities (we recall that V is a decreasing function in (0, oo)) 

°° °° 1 fks 1 poo 

VeffM = y~ v(ks) <JZ- y w dt = - / y w dt ' 



and 

1 r(k+l)s 



V eS (s) = J2 V(ks) > £ / y « dt =- y M 



dt. 
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and for a more detailed comparison between the internal stresses a> n { obtained from our derivation 
and the models proposed in the engineering literature. 

1.6. Related mathematical work on discrete-to-continuum transitions. The model of 
this paper lies halfway between one and two dimensions. Written as (1.1), it is a one-dimensional 
system, and an example of the general class of two-point interaction energies. There is a large body 
of research on this type of energy, which roughly falls into two categories. When the interaction 
energy is superlinear at infinity, the system models the behaviour of elastic solids, and examples of 
Gamma-convergence of such functionals are given in [12, Th. 1.22] (see also [1] and [2]). When the 
functional is bounded at infinity with a global minimum at finite distance, such as in the case of the 
Lennard- Jones potential r i->- r~ 12 — r 6 or the Blake-Zisserman potential r i-> min{r 2 , 1} [8], such 
two-point interaction energies lead to models of fracture (see e.g. [10], [9] and [13]). The functional 
V in (1.2) is neither of these, being purely repelling and convex away from the singularity. While 
the methods that we use are inspired by the general works in this area, we know of no work that 
deals specifically with this type of functional. 

There are various previous works that focus on the behaviour of minimizers rather than on the 
functional. The early work by Eshelby, Frank and Nabarro [18] mentioned before studies the case 
of a single row of dislocations (h = oo) and proves rigorously the asymptotic distribution of the 
dislocations. Hall [24] studies the wall setup, chooses the specific regime j3 n <~ n -1 / 2 , and proves 
convergence of stationary states using formal methods. Finally we should mention the numerical 
study [15] in which the correct asymptotic scaling of the regime 1/n <C f3 n <C 1 was already found. 
Mcsarovic and collaborators [6, 33] derive a continuum dislocation model from the discrete wall 
setup by means of a two-step upscaling: first the dislocations are smeared out in the slip plane 
and then in the vertical direction. Upscaling in the two directions separately, though, produces a 
significant error (referred to by the authors as "the coarsening error" ) that needs to be corrected 
by adding an ad hoc term to their continuum model. 

At the same time, the structure of the walls in Figure 1 is an attempt to make some progress in 
the problem of upscaling two-dimensional collections of dislocations. This is a hard problem, and 
the main difficulty can be recognized as follows. If we consider a field of edge dislocations in two 
dimensions at points {x^}^ C K 2 , and formulate the corresponding empirical measure on M 2 , 



The function V e d gc is singular at the origin, and therefore a simple weak convergence of /ijv in the 
sense of measures to some \i does not allow us to pass to the limit. 

To make things worse, c^V^jge takes both signs along the line x\ = constant. This indetermi- 
nacy causes a phenomenon of cancellation, and surprisingly this cancellation can be complete [36] : 
if we consider a continuous vertical line of smeared-out edge dislocations (i.e. the limit of a wall 
when /nO), then the total force exerted by this continuous wall on any other edge dislocation 
vanishes [36] . This cancellation is the reason why the tails of V decay exponentially, even though 
V dge only decays logarithmically. 

Because of the multiple signs of d Xl V e( ige and this cancellation, also a more advanced argument 
along the lines of [39] does not apply. Indeed, the results of this paper show how the relative 
spacing in horizontal and vertical directions has a major impact on the limiting energy. This 
relative spacing, the aspect ratio of the lattice of dislocations, is weakly characterised by /3„, 
which we discuss below. 




then the interaction energy for this system is essentially 




V r ed ge (x- y)/j, N (dx.)fi N (dy), where V edee ({x 1 ,x 2 )) 




1.7. Comments. In this section we collect a number of comments on the discrete model and on 
the results of this paper. 
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Conditions on V. While we perform the calculations in this paper for the exact functional V 
in (1.2), with minor changes the results can be generalized to any function V satisfying 

(1) V : K — > R is non-negative, even, and convex on (0, oo); 

(2) V has a logarithmic singularity at the origin; 

(3) V has exponential tails. 

On the choice of T- convergence. Our T-convergence result implies convergence of minimizers, 
and is stronger in a number of ways. For instance, T-convergence of E^ implies that E^ + F 
also T-converges whenever F is continuous. This allows us to deduce a similar convergence result, 
for instance, for a functional of the form 

n n / ~ ~ \ n 

1=1 3=0 V ' 1=1 

for any continuous function /, allowing us to consider more general, non-constant forcing terms. If 
in addition liminf x _ ) . (X) f(x) = +oo, then a similar compactness result also holds. Note, however, 
that such a functional obviously behaves differently under rescaling of the a;,. 

A second reason why T-convergence is a stronger result is the role that it plays in convergence 
of the corresponding evolutionary problems, i.e. the ordinary differential equations (1.13). That 
system is a gradient flow, and a method such as in [38] makes use of the T-convergence of £ (and 
other properties) to pass to the limit n — > oo in such a system. 

Connection between the limit junctionals. The transitions between the five different limiting 
functionals of Theorem 1.1 are continuous. For instance, if in E^ in (1.7) we take the limit c — > oo, 
then s i ^ cV(cs) converges to (J V)5, and we recognize the corresponding single integral in (1.8). 
In the case of E^\ in the limit c^Owe approximate V e s(s) by its leading order Taylor-Laurent 
development at the origin, which is (l/2s) j R V by (1.14), upon which E^ becomes equal to E<- 3 \ 

Similar transitions exist from E^ to E^> in the limit c — > 0, and from E^ to E^ in the limit 
c — > oo. 

Boundary layers. Figure 5 shows a good match over most of the domain, with a sharp boundary 
layer near the origin. The reason for this boundary layer can be recognized in the fact that 1 /nf3 n rts 
0.08 is about one order of magnitude smaller than the domain of the density. Such boundary layers 
are well known in the theory of interacting particles with next-to-nearest neighbours (see e.g. [9]), 
and we believe that the effect here is similar. Note that in Figure 6 the boundary layer is thinner, 
and indeed there l/nf3 n w 0.006. 

Generalisations. Baskaran et al. [6] and [33] study the same setup with arbitrary angle between 
the slip planes and the obstacle. They point out that orthogonal slip planes are a special case 
among all angles, and an obvious avenue of generalization is to understand the general case. 
Other generalizations include dislocations of multiple signs, creation and annihilation effects, and 
convergence of the evolution equations. 

1.8. Organisation of this paper. In Section 2 we prepare the stage for the main proofs, by 
introducing equivalent formulations for the rescaled energies and a characterization for the lower- 
semicontinuity of the limit functionals. Section 3 is devoted to the proofs of the five cases of 
Theorem 1.1, and Theorem 1.2 and Corollary 1.3 are proved in Section 4. 

2. Preliminaries 

In this section we collect a number of preliminary steps leading to the proof of Theorem 1.1. 
We start with rewriting the discrete functionals in a number of different, equivalent forms. In 
Section 2.3 we derive an equivalent characterization of the weak convergence of measures, and in 
Section 2.4 we characterize the lower semicontinuous envelope of functionals of the form J f(u'). 
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2.1. Notation. Here we list some symbols and abbreviations that are going to be used throughout 
the paper. 

f2 domain [0, oo) 

C one-dimensional Lebesgue measure restricted to fl 

A4(Q) non-negative Borel measures on of mass 1 

Cb(A) continuous and bounded functions in A C R 

IImIItv(A) total variation of fj, € M(£l) in A C 

dv/du Radon-Nikodym derivative of v with respect to fi 

fi„ n n fi„ ® fi n without the diagonal terms (see (2.1)) 

E^ k \ k € {1,2,3,4,5} discrete energies (see Theorem 1.1 and Section 2.2) 

E (k) , k € {1,2,3,4,5} limit energies (see Theorem 1.1) 

Also, we write e.g. J n f d/i instead of J °° / dfi, since the latter is ambiguous when /i has an 
atom at zero. 



2.2. Rewriting the functionals. The continuum limit functionals E^> and E^ are convolu- 
tion integrals, and this suggests reformulating the corresponding functionals at finite n also as 
convolution integrals. For given \i n = h Y^ii=i ^ t > we define the measure fi n M fj, n as the product 
measure fi n ® [i n without the diagonal: 

^ n n 

^„H M „(A) := -2 ^2^2S (XiiXj) (A) for any Borel set Ac n 2 . (2.1) 

n »=1 3 = 1 

Omitting the diagonal does not change the limiting behaviour: 
Lemma 2.1. If /j„ — jj,, then [i n M n n — ^ ji ® /U. 
Proof. Take ^(7 t (!l 2 ). Then 

/ tpd/J, n M [l n - I ipd^tg) fj, — tpd(/J, n M H n - /i„ ® H n ) + / <pd((jLn<8) fJ, n - fj,® fj,). 

Jn 2 Jn 2 Jn 2 Jn 2 

The second term on the right-hand side converges to zero since fj, n — ^ fi, and the first is bounded 
by H^Hoo/n and therefore also converges to zero. □ 

With this notation we can write in a number of different, equivalent forms: 
EW(x u . . .,x n ) = £(x u . ..,£„) + ^ (log27rn 2 ^ - 1) 

^ n n \ n 

= ^ E E Mn 2 (3t( Xi - Xj )) + - *i 

i—1 j—1 i—1 



= ^2 E E Vn(n 2 (3 2 n (x ]+k - Xj )) + - E ^ 

k=l j=l i=l 

= \ I V n {n 2 [3 2 n {x - y)) n n M ii n {dxdy) + [ Xf i n (dx). (2.2) 
2 Jf22 Jn 

Here V„(s) := K(s) + 7r~ 2 (log(27rn 2 /3 2 ) — 1) is a renormalized energy, obtained by removing a core 
energy from the energy density V. 
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Similarly we rewrite 

^ O^l > • ■ • ! x n) ^ {%1 > • • • i — ^ (x\ , • ■ ■ , X n ) 

Pn p i ~ ~ \ 

« n n— fc 
Pn 



XI H V(n0„(arj + k - Xj )) + -^Xj 
fe=i i=i j=i 

- / / V"(n/8 n (a; - y))Mn ^ Vn{dxdy) + / (da;), (2.3) 



and 



°V!r r n^ fe=l j=l \ \ / / J=1 

2.3. Convergence concepts and compactness. As already discussed in the introduction, there 
are two natural ways of describing the positions of a row of dislocation walls: 

(i) The position xf as a function of particle number i. One can make this formulation slightly 
more useful by reformulating it in terms of increasing functions £™ : [0, 1] — > fi, such that 
£ n (i/n) — xf, with linear interpolation. 

(ii) A measure fj, n = ± Yh=o S x»- 

In the introduction we mentioned the formulation in terms of measures as the basis for conver- 
gence results. However, in the proofs it will sometimes be useful to use the formulation in terms of 
functions £™ . Since we intend the resulting T-convergence to be independent of which formulation 
we choose, we choose a single concept of convergence and formulate this equivalently for £™ and 
for \i n . This is the content of the next theorem. 

Theorem 2.2. Let (xf) be a sequence of n-tuples such that Xq = and x" < xf +1 for every n 
and for every i = 0, . . . , n — 1. Let £™ : (0, 1) —¥ R + be the affine interpolations of x n , i.e. 

TO := + n(s? +1 /or S e(^,^), (2.5) 

and define the measures pL n e A4(Cl) by 

1 ™ 

M« : = - Yl s x?- (2-6) 
n 

i=i 

Then the following convergence concepts are equivalent: 

(i) £™ converges to £ in BV(0, 1 — 6) for each < S < 1 (we indicate this as 'convergence in 
BV loc (0,l)'); 

(ii) ji n converges weakly to \i. 

If the limit function £ is strictly increasing and a.e. approximately differentiable, then it is related 
to the limit measure /i by the formula 

" ( *» - mkrr (2 - 7) 

Finally, if 

1 " 

sup-V<<oo, (2.8) 
n n i= i 

i/ien i/ie sequences fi n and £™ are compact in this topology. 
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Note that the weak topology on the space of non-negative Borel measures A4(Cl) of unit mass 
is generated by a metric (see e.g. [4, Remark 5.1.1] or [7, p. 72]), and therefore there is no need to 
distinguish between compactness and sequential compactness. 

Proof. First we prove (i) => (ii). Let J" : [0,1] — > R + denote the piecewise constant function 
such that £ (s) = x™ for s e (^-, ^] , for every i = 1, . . . , n. We have for sufficiently large n, 

r is _ , r»(i-5)i , 

0</ (?(8)-C{s))ds<- Yl « + l-^)<-Hri|w(0,l-5/2) rl -^ O 0, 

« i=0 n 

so that J™ — > £ in £i oc (0, 1) and, after extracting a subsequence without changing notation, £™ — > £ 
pointwise a.e. 

Let now p g Cf,(]R) be a test function (for the weak convergence of measures); then 
<p(y) n n (dy) = -J2 - -V(0) + £ / ds = + / 



and since £„ — > £ a.e., 

By the uniqueness of this limit the whole sequence [i n converges. By defining \i g .M(Q) through 



lim f\(f(s))d S = [ <p(t(s))ds. 

n ^°°J Jo 

t the whole sequence fi n converges. B; 

/OO ,1 
<p(y)v(dy)= / <p(Z(8))da, (2.9) 
-oo JO 

we have proved that fi n — 1 /U. 

The identity (2.9) expresses the property that is the push-forward under £ of the Lebesgue 
measure ds on (0, 1). It follows by [4, Lemma 6.5.2] that whenever £ is strictly increasing and a.e. 
approximately differentiable, then 

" {dy) = wkrr 

Next we prove (ii) (i). Since \i n is assumed to be of the form (2.6), we can construct the 
positions x™ and the linear interpolation £™ as above. The convergence of \x n implies that the 
sequence \i n is tight, which implies in turn that for each 6 > 0, 

sup sup x™ < oo, 

n i:i/ n <l-S 

and therefore that sup„ £™(1 — S) =: M < oo. 

Therefore, since £"(0) = by (2.5), we have the bound 

,1-5 ,1-5 



/" (l-s)(£ n )'(s)ds = S£"(l-(5) + / C(s)ds<M. (2.10) 
o Jo 



Therefore, using the monotonicity of x n we have that 

M> f 1 \l-s)(C)'(s)ds>5 f S \(C)'(s)\ds. 
Jo Jo 

This provides a uniform bound for (£")' in £ 1 (0, 1 — (5), and by integration also a uniform bound 

on £ ra in L 1 (0, 1 — 8). Hence the sequence (£") is equibounded in VF 1,1 (0, 1 — 5), and therefore 

converges in L 1 (0, 1 — 5) and weakly-* in BV(0, 1 — <5) to a function £ g BV(0, 1 — 5). 

Finally, the compactness of the sequence [i n follows from the tightness implied by (2.8) and the 

estimate 



| a; | jJL n {dx) 



r 1 i ™ 

/ d* <-£*?■ 

J( > n i=l 



□ 



Remark 2.3. Note that the limit function £ introduced in the previous theorem is increasing, since 
it is the pointwise limit of a sequence of increasing functions. 
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2.4. Lower semicontinuity and relaxation. This section is devoted to a lower semicontinuity 
result for functionals denned on the space of special functions with bounded variation. More 
precisely, the next theorem provides an integral representation for the relaxed functional in a 
special case. 

Theorem 2.4. Let f : (0, oo) R be a convex and decreasing function such that lim t _ ) . 00 f(t) = 0. 
Let F : BVi oc (0, 1) — > R U {oo} be the functional defined as 

. = I J f{u') dt if u eW 1 ' 1 (0,1), u increasing, {2 11) 

y+co otherwise. 

LetH denote the lower semicontinuous envelope of F (relaxation of F) on BVi oc {0, 1) with respect 
to the BV\ oc {0, 1) -convergence defined in Theorem 2.2. 

We introduce the functional T : BVi oc (0, 1) — > M defined, for u increasing, as 

F{u) := / f{u')dt, (2.12) 
Jo 

where v! denotes the absolutely continuous part (with respect to the one- dimensional Lebesgue 
measure) of the measure Du, which is the distributional gradient of u. Then we have 

T=ri. 

Proof. We first note that by construction T < F on BVi oc (0, 1). Since T is lower semicontinuous 
with respect to strong convergence in Ll oc (0, 1), by e.g. [3, Proposition 5.1-Theorem 5.2], it follows 
that 

T<n. 

For the opposite inequality we need to show that, for a given u e BVi oc (0, 1), u increasing, there 
exists an approximating sequence (u e ) C VK 1,1 (0, 1), u l increasing for every £, such that u e — > u 
in L\ oc and 

limsup f /((?/)') dt< [ f(u')dt. (2.13) 

£-s-co Jo JO 

For the construction of the sequence [u l ) we proceed as follows. We first approximate the 
distributional gradient Du of u with L 1 functions, say w e , with respect to the weak convergence 
in measure. Then we construct approximations u l as (properly defined) anti-derivatives of w l 
and will be therefore in W 1 ' 1 by construction. This argument is strictly one-dimensional, since it 
makes use of the property that every function is a gradient. 

We now go through the details of the proof. 

Step 1: Approximation of Du with L 1 functions. We decompose the distributional gradient 
Du as Du = u' + D s u into its absolutely continuous part and singular part with respect to the 
Lebesgue measure. Since u is increasing, both u' and D s u are non-negative measures (being 
mutually singular). We notice that the absolutely continuous gradient u' (identified with its 
density with respect the Lebesgue measure) is by definition a nonnegative L 1 -function; therefore 
it is sufficient to approximate the singular measure D s u with nonnegative functions in L 1 . Let 
(g e ), with g l <G i 1 (0, 1) be such an approximation and define 

w l :=u'+g l ; (2.14) 

then w l e L 1 (0, 1), w l > a.e. and w l — Du weakly in measure. 

Step 2: Approximation ofu. We notice that, for the construction of the approximating sequence, 
we can assume that S u = 0. Indeed, let us assume instead that S u ^ 0; by the locality of the 
argument we are going to use, it is not restrictive to assume that S u = {t*}. 

We define continuous approximations of u as 

'u(t) if t e (0,t* -e)U (t* +e, 1) 



u e (t) :-- 



u{ t * - £ )+ u(t'+e)-u(e-e) ^ _ t * +e) if t g (j* _ £; t * + £ )_ 
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Then, clearly, 

(u'(t) if te (0,r-£)U(i*+e,l) 
«*(*) == { 

( u{t +£) -" (t " £) if te (t*-e,t* +e). 
For the approximating sequence w £ we have 

fl /• /•**+£ 

2e 



f(u' s (s))ds = I f(u'(s))ds + I f ( " (r+£) " (t * £) ) rf, 



J(0,l)\(t*-e,t*+e) 

f(u'(s))ds + 2sf 



(0,l)\(t* -e,t*+e) 

< [ f(u'(s))ds + 2ef 
Jo 



u(t* + e) - u(t* - e) 



2s 



u(t* + e) - u(t* - e) 



2e 



The decay at infinity of / implies therefore that 



limsup / f(u' E (s))ds < / f(u'(s))ds. 

6^0 Jo Jo 



Therefore we can assume that u is continuous. 

We define the primitive of the function w l defined in (2.14) as 

u\t) := u{0) + f w e {s)ds. 



It follows that u e £ T4 /1,1 (0, 1), u e (0) = u(0) and (u e )' = w e , which converges weakly to Du in 
measure. 

Since (u l ) is bounded in W 1 ' 1 , then it converges weakly in BV to a function v G BV(0, 1). By 
the weak convergence of (u e )' in measure it follows that Du = Dv and therefore, since v(0) = u(0), 
that u = v. Hence, we have constructed a sequence (u e ) C 14 7l ' 1 (0, 1) such that u l — > u in BT-^O, 1), 
and hence in £>Vj oc (0, 1). 

Step 3: Upper bound for the energies. Since (u e )' = u' + g l and g e > we have by construction 
that 

f < /' /(«')<*« 

Jo Jo 

for every £, since / is a decreasing function. The bound (2.13) follows immediately. □ 

3. Proof of Theorem 1.1 
We separate Theorem 1.1 into the five different cases, and state and prove each case separately. 

Theorem 3.1 (Case 2, first critical regime: I3 n <~ l/n). Let c n := nj3 n — > c > as n — > oo. -For 
£/iis case the functional E^ in (2.3), which can be rewritten as 



E nH^n) = y V(c n (x - y)) fj, n Kl fj, n (dxdy) + J x [i n (dx 



T-converges with respect to the weak convergence in measure to the functional defined for 
/i E M(fl) as 

£ (2) (M) := % [[ V(c(x - y)) n{dx)^{dy) + [ x^dx). (3.1) 
^ J Jn 2 Jn 

In addition, if En \^n) is bounded, then \i n is weakly compact. 

Proof. The compactness statement is a direct consequence of Theorem 2.2, since the interaction 
potential V is non-negative. The remainder of the theorem we first prove under the assumption 
that c n = 1. 
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Liminf inequality. Let \i e A4(fl) and let fi n be a sequence of measures of the form \i n — 
i Y^i=i sucn that fi n — ^ weakly in measure. Since V > is lower semicontinuous on R 2 and 
fi n M ji n — ^ p, ® ji by Lemma 2.1, 

liminf i// V{x - y) p n M p n {dxdy) >\ ( ( V (x - y) fj,(dx) fj,(dy) . 

For the second term we have a similar bound, and therefore 

hminf4 2 )( M „)>£;( 2 )( M ). 

Limsup inequality. It is sufficient to prove the limsup inequality only for a dense class, 

A := j/x e .M(fi) : supp,u bounded, /i < £, and ^ e 

This set is dense in A4(fi), and for any /i e M(Q) with E^ 2 \p) < oo an approximating sequence 
(/ife) C A can be found such that /Ufe — /i and E^ 2 \/j,k) — > E^ 2 \p). This can be seen, for instance, 
by defining 

Hk(dx) = pk{x) dx with pfc(x) = fc/n([x, x + 1/fc)). 

Then by Fubini, 

/•/■ /•/■ rX+l/k ry+l/k 

/ / V(x - y) n k (dx)fi k (dy) = k 2 / / V(a; - y) / n(d£) / /x(djy) dxdy 
J J a, 2 J Jn 2 Jx J y 

= [[ [ f V k 2 V(x-y)dydxp(d£)p(dr 1 ). 

JjQ 2 J{i-l/k) + J(r,-l/k) + 

Now one recognizes in the inner two integrals the convolution of the function 

(x, y) V(x - y)xn{x)xn{y) 
with the characteristic function of the square [0, 1/fc) 2 , so that the expression above converges to 

V(x - y) n(dx)n{dy). 

o 2 

This shows that it is sufficient to prove the limsup inequality for all /i € A. 

Take such a measure fi € A with Lebesgue density p G L°°(Q), and construct an approximation 
M« = Tl Z)"=i <W b y defining the points Xj by 

/ p(x) dx = — . 
Jo n 

Then 

Since supp p is bounded, all x™ are uniformly bounded, and 

lim / x n n (dx) = \ x p(dx). 
™^°° Jo Jo 

Turning to the convolution term, for fixed m > we write 

\ jj^Vp n Mp n = X - ff n2 ( VAm ) Vn^Pn + \ JJ^(V - (V Am)) fj, n ^ /j, n . 

In the first term the function V A m is bounded and continuous, and this term therefore converges 
to 



\ JJ (VAm)ii®ii< l - JJ Vn®n. 
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If X m > solves V(X m ) — m, then we estimate the second term by 

Iff {V ~{V Am))fJL n Mll n < \ [[ VHn^Vn 
1 J JQ? Z J J {\x-y\<X m } 

^ n n—k 

-2 E E v ^+k - x ?) 1 {\'? +h -^\<x m } 



k=l 3 = 1 

L^m«|M|ooJ n—k 



^ h E E - »") 
fc=i j=i 

< - V v(- 

LX m n||p|Uj 



E n-)* 

7 i ^ - n _ m 



DC 



fe=l 

V(s)ds. 



o 



Therefore 

limsup£;i 2 )( M „) < £< 2 )( M ) + HplU / V(s)ds. 

n—>oo Jo 

Since m > is arbitrary, and since lim m _ ) . 00 X m = 0, this proves the limsup estimate 

hmsup4 2 >( M „)<£( 2 >( M ). (3.2) 

n— >oo 

In order to allow for c n ^ 1, we define the scaled measure 
with which 



n <= i 



E n } (Vn) = % / / ^(x - H (i n (dxdy) + — [ xfj, n (dx). 
* J Jn 2 °n Jn 

The two prefactors in this expression do not change the arguments above, and upon back- 
transformation the result of the theorem is found. □ 

Theorem 3.2 (Case 1, subcritical regime: /3 n -C 1/n). Let f3 n > be a sequence such that 
n(3 n — > asji^oo. Then the junctionals En defined in (2.2) Y -converge to the functional E^ 
defined on measures [i G A4 (O) as 

^ (1) (M):=-^ JJ log\x-y\ii(dy)iJi(dx) + J xiM{dx). (3.3) 
In addition, if En\[i n ) « s bounded, then [i n is weakly compact. 

Proof. Compactness for the measures This is the only one of the five cases in which the 
compactness is non-trivial, since V n takes both signs; therefore a bound on E n ^ does not translate 
directly into a bound on the second term — ^2 x i- However, by combining the first two terms, such 
a bound can be obtained, as we now show. 
First we show that 

V(t) > V(t) := 1 ~ 1 ° S 9 27r|t| for all t ± 0. (3.4) 
w 

This follows by remarking that for t > 

V'(t)-V'(t) = -^- + ±->0, (3.5) 
sinh 7rf tt t 
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and using the expression 



V(t) = 



2t 



ir(e 27Tt - 1) 7r 



-^log(l 



-2lTt\ 



we compute that for t > 



]imV(t) - V(t) = lim 
*4.o tio 



^ T y + i{- 1 ^ 1 - e ^- 1 + 1 ^ 2 ^}] 



o. 



(3.6) 



From (3.5) and (3.6) we deduce (3.4). 

Therefore the renormalised interaction energy V n satisfies 



V n (n 2 p 2 n t) = V(n 2 p 2 n t) + 



log(2™ 2 /3 2 ) - i 



7T^ 



> 



l-log(27m 2 /3 2 |i|) log(27m 2 /3 2 ) - 1 



+ 



= 2 log 1*1 



(3.7) 



Note that for alii ^ 

V n {n 2 Plt) + \\t\>\\t\'l~ 2 log |t| > ± (l - log A) > o. 

Let (U„ be a sequence of measures of the form [i n = ^ ^z™ such that En\fj, n ) is bounded. 

We now estimate 



^ 1) (Mn)-^< 



2—1 j — 1 i— 1 j — 1 



1 " 
+ 2^E ; 



i=i 



1 n 1 

^EZ [^(n 2 ^(»i ~ ^)) + 2^ - ' 
i=i j=i 

1 ™ 



1 

2n 



E** 



i=i 



The boundedness of E n (fi n ) and Theorem 2.2 then provide compactness of the sequence fj, n . 

Liminf Inequality. Let now /j, n be a sequence of measures of the form [i n = ^ Y^i=i that 
converges weakly to it, and note that by Lemma 2.1, fi n M fj, n — >■ /x ® /i. By (3.7) we have the 
bound 



E^{fin) >-^JJ^Jog\x-y\fi n ^ (i n (dxdy) + J x fj, n (d. 



'x) 



S So, 2 



1 1 1 If 

—rlog\x-y\ + -(x + y) {i n M ii n (dxdy) + - I x^ n (dx). 
2ir z 2 J n J n 



(3i 



The function between brackets is lower semicontinuous, and by a similar argument as we used for 
the compactness above it is also bounded from below. Therefore the right-hand side in (3.8) is 
lower semicontinuous with respect to weak measure convergence, and therefore 

liminf E^\fi n ) > --^ [f log \x - y\ /j,(dy)n(dx) + [ x n(dx) 
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Limsup inequality. For the construction of a recovery sequence we first prove a second 
inequality on V for t > 0: 

7r(e 27rt - 1) 7T 2 



-2-rtt 



< 



1 



< 



7T Z 
1 



-2nt + log(e 27It - 1) 



1 + 2irt - log 2irt 



from which follows the estimate for all i 7^ 0, 



V n (n 2 p 2 n t) = V(n 2 p 2 t) + 



log(27m 2 /3 2 ) - 1 < 1 



: (2irn 2 f3 2 n -\og\t\). 



(3.9) 



7T^ IT 

The remainder of the argument follows largely the proof of Theorem 3.1. Given a similar limit 
measure [i and approximating sequence p n , we estimate 



< 



V {n 2 P^(x - y)) p n M p: n (dxdy) + / x fi n (dx) 

Jn 

Iff f n 2 p 2 
2~^JJ l0 S\ x ^y\l^n^Hn(dxdy)+ J x /j. n (dx) + -. 



Decomposing — log |a; — y\ into a part that is bounded and a remainder, as in the proof of Theo- 
rem 3.1, and repeating the corresponding estimate, one can show that the right-hand side converges 
to E^(fi). This proves 

limsup E^(ti n )<E^ (12). 

n— >oo 

□ 

Theorem 3.3 (Case 3, intermediate regime: 1/n <C p n <C !)• Let f3 n > be a sequence such that 
P n — > and nP n — > 00 ; as n — > 00. T/ien i/ie Junctionals E n 3 ^ defined in (2.3) V -converge with 
respect to weak measure convergence to the functional E^ defined on measures (1 e A4(fl) as 

^ V I /" p 2 (x)dx+ / xp{x)dx if jji = pdx, 



+00 



otherwise, 



which is the same as 



V 



-ds - 



£(s)ds, 



(3.10) 



(3.11) 



/o £'(*) 

w/ien written in terms o/£ £ -BV/ oc (0, 1), £ increasing, and /1 and^ are linked by (2.7). In addition, 
if En\nn) ^ bounded, then fi n is weakly compact. 



Proof. Again the compactness statement follows from Theorem 2.2. 

For the liminf inequality we will make use of the expression (2.3) for the energy, i.e., 

1 



E n 3 \ 



■ nP„ 



Q 2 



V(nP n (x - y)) ii n M jj, n {dxdy) + / xfi n (dx). 



We will prove that for any sequence \i n — 1 (j,, 

\immiE^(p n )>E (3 Hfi). 



(3.12) 
(3.13) 



Take a sequence \i n — >• (i such that E^\p, n ) remains bounded. Since the second term of 
En (Mrt) is bounded, Theorem 2.2 guarantees that there exists a measure fi such that [i n — »■ /1 
in measure, at least along a subsequence, and we switch to that subsequence without changing 
notation. The support of the limit measure fi lies in O = [0, 00) by the definition of the extended 
measures \i n . We split the rest of the proof into three steps. 
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Step 1: Rewriting the energy in terms of convolutions. Let V n (t) := nf3 n V(n/3 n t); we claim that 
V n converges to (J R V)5q in distributions. Indeed, let ip e C^°(K); then 



lim / V n (s)^(s)ds = lim / V(t)i/> ( ^~ 

™^°°JR n ^°° JR \nPn 

which proves the claim. 



dt=[ V(t)dt V(0) 



Now we use the fact that V = W * W, where W = U and U = Vv, as proved in the Appendix 
(Subsection A.l). In addition to W we will also use its truncation W m :— min{ W,m} for any 
fixed m > 0. 

Let W n be defined as W n :— U n , with U n = \fVn- By the scaling properties of the Fourier 
transform it follows that W n (t) = nj3 n W{nj3 n t). Similarly we define W™(t) := n/3 n W m (n(3 n t). 
We note that, like in the case of V n , the distributional limits of W n and W™ are (J R W) do and 
(J R W m )S . Moreover, 



f w{t) dt = w(o) = Vy(o) = J / v{t) dt. (3.14) 

JR 



Note that since V = W * W, 

V n {x-y)= / W n (z-x)W n (z-y)dz, 



and therefore 



2 / / K(z - y)/i„ Kl /j, n (dxdy) = — ^ K(^+fe - Zj) 
^•^° 2 n fc=i j=i 

= ~EE / W n (z - x j+k )W n (z - xj) dz. 
71 



k=l j=l 



We then estimate 



\ [ (C*n) 2 4EE/ ^(«-^ + fc)W n ro ^-^)^+iE / W™(z- Xj fdz 

JR. fe=1 J=1 JR J = 1 J 

^-EE / w„(^-^ +fc )w„(^-^)^+-nw™iii 

= lf[ 2 V n(x- </K B Mdardy) + y IIWHI2 (3-15) 
< C. 

Therefore we obtain weak convergence in L 2 (M) along a subsequence of W™ * /u„ to some 

feL 2 (R). 

Step 2: Identification of f . In order to find the relation between / and fi we compute the 
distributional limit of the sequence * [i n . Let tf> € Cq°(IR) be a test function; then we have 

lim f {W™*ii n ){x)i>{x)dx= lim / {W? * ^){x)^ n {dx). (3.16) 

Note that W™ * V "— ? (/ R W m ) ip uniformly, since W™ * ip converges to (/ W m ) ip strongly in 
ff^R). The uniform convergence of W™ * ip, together with the weak convergence of fi n to /i, 
guarantee that the limit in (3.16) is: 

lim / (W™ * /j, n ){x)ip(x) dx = lim / (W™ * ip)(x)n n {dx) = ( [ W m ) f ip{x)n{dx). (3.17) 

JR n ^°° JR \Jr ) JR 
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Therefore, by the uniqueness of the limit of (W™ * p, n ) we deduce that / = (J R W m ) p. Hence, p 
is absolutely continuous with respect to the Lebesgue measure with a density in L 2 (M), i.e., there 
exists p e L 2 (R), p > a.e., such that p = pdx. 

Step 3: Lower bound. From (3.15) it follows that 

liminf Iff V n (x- y)p n m p n (dxdy) >\( f W m \ f p 2 (x) dx. 
Taking on both sides the supremum over m > we find 

liminf Iff V n (x- y)p n ® p n {dxdy) >\( f w] f p 2 (x) dx 

= l([ V )[ P 2 (x),dx, (3.18) 

where in the last equality we used the relation (3.14). For the second term of the energy we have 
that, for every M > 0, 

r- r- pM 

liminf / x p n (dx) > liminf / x p n (dx) — f xp(x)dx, (3.19) 

™^°° JR J[0,M] JO 

so that, taking the supremum on all M > 

liminf / x p n (dx) > I xp(x)dx. (3.20) 
n-J-oo J M J n 

In conclusion, from (3.19) and (3.20) follows the inequality 

liminf E^\p, n ) > \ ( f V] f p 2 (x)dx+ f x p(x) dx = E^(p) , 
2 \J m J J n J n 

which is (3.13). 

We now continue with the proof of the limsup inequality: for each £ <G -BV[ oc (0, 1), there 
exists a sequence (x n ) n of n- vectors (x™, . . . , x") such that 

limsup EW{x?,...,x?l) < E^(0 = l(J R V) J -^ds + J Z(s)ds. (3.21) 

By Theorem 2.4 we can assume without loss of generality that £ € W ' (0, 1). 

So, let £ e VK 1:1 (0, 1) be an increasing function such that i?( 3 )(£) < oo. We can assume that 
there exists e > such that £' > e uniformly on (0, 1). Indeed, we can otherwise approximate £ 
by the sequence £ e (i) := £(t) + et. Clearly £ £ — > £ in T'F 1,1 as e — > 0; moreover for the absolutely 
continuous part of the distributional gradient of £ £ we have that = + £■ Hence ^ > e, since 
£ is increasing. Also, since > £' and f(t) — 1/t is decreasing we have 



I W) dt ~l W) dt > 



so that 

limsupS (3) (4) < E {3) (0- 
Therefore, from now on we can assume that £' > e for some e > 0. 

For every n & N we define the piecewise affine function £™ and the points x™ by £™ ( A ) := a;™ := 
£ («)■ Clearly the sequence £" converges to £ strongly in W 1 ' 1 . We consider the energy for this 
sequence, 

P n n—k n 

^i 3) w, • • • , o = £ E E ^ - <)) + „ E <• 

fc=l i=0 i=0 
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As the second term of the functional is the Riemann sum of the integral of £™ in (0, 1), we focus on 
the first term. By the convexity of the energy density V we have that, using Jensen's inequality, 

( i + k \ i+k 

for every k = 1, . . . , n and every i = 0, . . . , n — k. Since 

? I L " v (/?nA:(r) ' (s)) ds - f v (Mn'W) ds for evcr y * = 1. ■ ■ ■ , «. 

we have the following estimate for the first term of the energy: 

q n n—k n „i 

- E E y (m^* - *?)) <^E / ^ y w) & 

" fc=l *=o fe^- 70 

= jf (e^R^" r )/(s) S ^ (*0W (*))) d*. (3.22) 

We define <5 ra (s) := fi n (£ n )'( s )- Since by assumption (£")' > e a.e. in (0,1), it follows that 
nS n (s) — > oo for a.e. s € (0, 1). Moreover, since £' is finite for a.e. s, S n (s) — » 0. It follows that for 
a.e. s, the expression 

n 
fc=l 

is a Riemann sum for the integral f °° V(t)dt = (1/2) J R V. Therefore letting n — > oo and by virtue 
of (3.22), we have the following bound for the energies ^ ; 

limsu Pj B( 3 )K,...,x^)< if/ v) [ j^ds+ [ £{s) ds = {Q. 
This proves (3.21). □ 

Theorem 3.4 (Case 4, second critical regime: n <~ 1). Le£ P n > be a sequence such that 
P n — > c > as n ^ oo. Then, as n — > oo, £/ie junctionals Eyf 1 defined in (2.3) T -converge to the 
functional defined in terms of measures /z e .M(fi) 6y 



^ ;= . ^(^y)p(-) d* + ^ *P(«) d* i/ M = Pdr, (3 23) 

[+oo otherwise, 

or in terms of increasing functions £ e £?V; OC (0, 1) as 

£ ( %) = c [ V eS (c?(s))ds+ [ Z(s)ds, 
Jo Jo 

where V e s(t) := X^fcli ,/ or every t G E. and fi and £ are linked by (2.7). In addition, if 

E^\fJ, n ) is bounded, then [i n is weakly compact. 

Proof. Again the compactness follows from Theorem 2.2. We first prove the theorem under the 
assumption that c = 1, and comment on the general case at the end. 

Liminf inequality. We will show that for every sequence (x™, . . . , x™) n of n-tuples, converging 
to £ in BV\ oc in the sense of Theorem 2.2, 

liminf 4 4 )(^..., <)>£(%). (3.24) 
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Take such a sequence (a;™, . . . , x™) n . We first rewrite the functional E^ in a more convenient way. 
For every fceNwe define the function V k (t) := V(kt). Hence from (2.3) we have 

fc=l j=0 \ n J i=0 

The expression in the argument of V k resembles a gradient, and we make the gradient term appear 
explicitly using the affine interpolant £" of the (see (2.5)). For fixed k € {0, . . . , n — 1} and 
£ € {0, . . . , n — 1} and for i<£<k+i— 1, we have 

n n -i 2— ^+fe— 1 n n -i i— ^+fe— 1 

*i+fc-*i_J £ X i+m+l ~ Xl+m = 1 V- (r) //- s+ ^^| (3 26) 

— /c \ ti / 



k y, 

n m=i-l n " m =i-t 



for every s € ( ^ , ) . Then 

i+fe-1 /-, i-^+fc-1 



n \ -. i-t-K— i /■i^- / i i— t-t-re— j. \ 



lyfe [ x »+fc x 
n I 

3=0 n \ m=-j 

Therefore, we can rewrite the first term in (3.25) in terms of the function x, as 

n ~ k 1 ( „n „n\ -i n-fefc-1 / 1 fc-l-j 



ds 



i=0 j=0 n y TO=-j 

fe-1 + 1 / fe-l-j 

lE/ ■ v i E(rr = 



1=0 

and the first term of the functional E n becomes 



fe-l I fe-l-j 



(is 



cis, 



i E E v P±V^ - E iE/ W E «">' • + 1) <- ^ 



fe = l i=0 \ n / fe=l j=0 



Now fix (5 > and note that by Theorem 2.2, £ n converges to £ weakly in SF(0, 1 — <5). We 
claim that for every fixed integer AT the following lower bound is satisfied: 



lim inf i#) (*?,... ,<)> f 5 V e ^'(s))ds+ f \{s)ds, (3.28) 

where the energy density Vjg is defined as Vjg (t) := J2k=i V(kt) for every t £ R. This claim 
implies the lower bound (3.24) by the arbitrariness of N and of 5. 

As in the proofs of earlier theorems we focus on the first term of the discrete energy E^; 
indeed, since £™ — > £ in L 1 (0, 1 — 5) (by the SF-convergence) , the bound on the second term of 
the energy in terms of the integral of £ on (0, 1 — 5) follows. 

Let N be fixed (independent of n). Then, for n> N, 



1 " ™- fe / ~n T n\ W fe-1 / fe-l-j \ 



fe=l i=0 \ n / fe=l j=0 



26 



DISLOCATION DENSITY 



We note that, since j and k run through a finite set independent of n, for every 77 > there exists 
an integer 1/(77) such that, if n > v(r]), then 

/ 1 _ -. — j 1 1 \ 1 / ' / 1 — jyi 



s H ) I gLs. 

n - 



for every j and for every k. Moreover, for every k and j, also the convex combination 

«,M.-iE~V(. + 2) 

m=-j 

converges to £ weakly in BV(r], 1 — 5 — rj), since it is bounded in W 1 ' 1 ^, 1 — 5 — 77) and has the 
same L 1 -limit as the sequence £™. Therefore, for every k and j 

liminf /' * \ k ((Zl j )'(s))ds> f 5 \ k {i'{ S ))ds, (3.29) 

since the integral functional in (3.29) is lower semicontinuous with respect to the weak convergence 
in BV, by e.g. [3, Proposition 5.1-Theorem 5.2]. In conclusion, 

JV fc-l r l-S-n r l-8 



fc=l j=0 ,/ '' 1/0 

= £/ V k ^( S ))ds+ d( S )d S 

k=1 Jr, JO 

= f 5 V V c %(£(s))d S + f \{s)ds, 

Jri JO 



and the claim (3.28) follows by the arbitrariness of rj. 

Limsup inequality. By Theorem 2.4 we can reduce to proving the existence of a recovery 
sequence for a function £ <G M /1 ' 1 (0, 1). 

Therefore, let £ <G 14 7l ' 1 (0, 1) be an increasing function such that E^^) < 00. For every neN 
we define the piecewise affine function £™ and the points xf by £ n (^) := xf := £ (^). The 
sequence £™ converges to £ strongly in W 1 ' 1 , and therefore also in BV\ oc . 

As in (3.25) we write 

1 n n-k / n _ ~n\ i n 
fe=l j=l \ n / j=0 

Since the second term of the functional converges to the integral of £ in (0, 1), we focus on the 
first term. As in the proof of the previous theorem the convexity of the function V k implies, by 
Jensen's inequality, 

n n — k / n _ „ra \ -, n n—k / \ n n—k „i±± 

fe=l j = l \ 71 / fe=l j=0 \ J n / fc=l J=0 J n 

(3.30) 

Since 

ra-fc , ^ ,.1 

V- / " V k (£'(s))ds< / V fc (£'(s))ds for every k = l,...,n, 
we have the following estimate for the energy: 

1 Tl n— k / n n \ n f\ (-1 

n J2E vk <E/ v k (?(s))d s < vMtWds, 

U k=lj = l \ n J k=l J ° Jo 
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which proves the desired inequality, 

limsup£;( 4 )K,...,<)<^ 4 )(C). 

General c. The case of general c — lim n _). 00 (3 n follows by rescaling, as in the proof of Theo- 
rem 3.1. In terms of the scaled measure 

1 ™ 

^ ■ i 

the functional E$f* reads 

E n i] (Pn) = ^ // V(ar- y)£„ H ji„(<to<%) + / a; 

* J JQ 2 Pn Jn 

Since /?„ is bounded away from zero and infinity, the same arguments apply, and we find that the 
right-hand side T-converges to 



Vcfi { —, z \ p(x) dx H — / xp{x)dx ifp^pdx, 
\p(x)J C In 



[+oo otherwise, 
where p and p are linked by 

/ ip(x) p(dx) = / tp(cx) p(dx) for all ip € Cb(f2). 

Back-transformation gives the T-convergence of the unsealed to the #( 4 ) defined in (3.23). □ 

Theorem 3.5 (Case 5, supercritical regime: !3 n — > oo). Lef /3„ > be a sequence such that 
I3 n — > oo as n — > oo. TTien i/ie junctionals E n ^ defined in (2.4) T -converge with respect to the 
strong convergence in Lj oc to the functional E^ defined for £ € BVi oc (0, I), £ increasing, as: 

£ (5) (0 = 11 ^ S)dS J/? ' (S) - 1 ^ °' e ' ° < S < (3-31) 
[ +oo otherwise, 

or equivalently, in terms of measures p linked to £ by (2.7), 

Eib)[p)= [^xp{dx) ifp<C, 

y +oo otherwise. 

Note that the inequality p < £ is intended in the sense of measures, i.e. p(A) < C(A) for all 
A C fi measurable. Equivalently, one can require that ^ <C C and d^i/d£ < 1. 

Proof. First of all, we define the sequence a n :— 5^1°g(f/^) and rewrite the energy (2.4) in 
terms of the new sequence, as 

E n 5 \ Xl ,...,x n ) = -Y J Y, V ^ na n( x j+k-x j )) + -Y J x r (3-32) 

/ nan k =ij=i 71 j= i 

Notice that li 2 n = | e 2Ta " and that a n > 1, since /3„ > 1. 

Liminf inequality. Let (x™, . . . , x") be a sequence of n- vectors such that the piecewise affine 
interpolation £™, as defined in Theorem 2.2 converges in BV\ oc (Q, 1) to some £. As for the other 
cases, the second term in the discrete energy E n 5 ^ is lower semi-continuous with respect to this 
convergence, and therefore we focus on the first term. 

The energy satisfies the trivial estimate 

n n—k n—1 
k=l i=0 i=0 
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since every term in the sum is nonnegative. Now, following the proof of the liminf inequality in 
Theorem 3.4 we can write 



n 

and therefore 



-V (a„ i+ \ j j - J i V(a n (C)'(s))ds, 

-, n-l „1 

- E V (na n (x? +1 - x?)) = / V(a n (CY(s))d S . 



i=0 

Hence, we have the following bound for the first term of the energy: 

D 27ra„ n n ~ k 



7T e 



2 na„ 

fc=l i=o 



EX>(™"(*?+* "*?))> / V n ((C)'(s))ds, 
u-1 .-n JO 



where V^(t) := f e a V(a n t). We claim that for any t n — > i, 

liniinfy rl (i„)>T/ 00 (t):={+ ° ^f^ 1 ' (3.33) 

10 if £ > 1. 

By e.g. [11, Prop. 2.2], this inequality implies that 

liminf^K,...,^)^^ 5 )^). 

n— ^oo 

To prove (3.33), we only need to show that if < t < 1, then V^(t„) — > oo. This can be easily 
proved in the following way. Since < t < 1 can be rewritten as t := 1 — r\ (for < 77 < 1), then 
for a sequence t n converging to t, we can assume that \t n — t\ < |, and in particular t n < 1 — |. 
Since is decreasing for every n, we have the bound V n (t n ) > V n (l — ^) for n large enough. 
Therefore 

liminf > lim V n ( 1 - - ) = 00, 

71— >oo n— >oo \ Z / 

which concludes the proof. 

Limsup inequality. We can once more invoke Theorem 2.4, since f(t) = (sc~Voo)(i) satisfies 
the assumptions of the theorem, and construct a recovery sequence only for increasing functions 
£ G W 1 ' 1 ^, 1) such that £( 5 )(£) < 00. 

By density we can further reduce to £ piecewise affine. 

Let us first assume that £ is linear, i.e., £(s) = (1 + £)s, for £ > (since i?( 5 )(£) < 00). 

The recovery sequence xf can be constructed in the following way. Let 5 n > be a sequence 
such that <5„ — > l\ then we set x™ := (1 + S n )^ for i = 0, . . . , n. The sequence (xf) is increasing, 
Xq = for every n, and 

E ^ (x?, . . . , *») =| e — E E ^ («*»(^* - *?)) + - E *? 

n fc=l «=o i=0 

27TO n n_ ^ /-I £ \ n 

= \ e —Y.Y.V Ml + ft,)) + ^ E *■ (3-34) 

™ fe=l i=0 i=0 

We claim that the sequence 5 n can be chosen so that 

„27ra„ " "- fe 



lim - — - {a n k{l + S n )) = 0. (3.35) 



fc=l i=0 



We focus on the term k = 1 in the sum in (3.35); we note that it is an upper bound for the other 
terms in the sum, corresponding to k > 2. Using the form of V for large values of t recalled above, 
we can rewrite the term k = 1 as 
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n-1 



e 



V V (a n (l + 5 n )) «K(1 + S n ))e- 2 ^ 1+5 ^ 



= (1 + <5„)e- 2 ™"S 

up to a remainder going to zero exponentially fast as n — > oo. Clearly, the sequence 5 n can always 
be chosen so that the last expression converges to zero (indeed, even if £ = we can choose S n — > 
such that a n 5 n — > oo). 

Therefore, the claim (3.35) follows directly, since every other term in the sum is estimated by 
the term k = 1; in conclusion, 

limsup^ 5 )^™) < -{! + £) = / x{s)ds = E^(x). 
To illustrate the general case of x piecewise affine we can reduce to x of the form 

x(t) = \ (3.36) 

[(i + e 2 )(t-t*) + (i + e 1 )t* ift>t*, 

where, for example < t\ < £2, and t* G (0, 1). Then assuming that nt* e N, the approximating 
sequence is defined as 

x n = \ ~ ~ (3 37) 

\{l + I2) - t*) + (I + li)t* iint*<i<n. 

Clearly x™ converges to x in L 1 . 
We claim that for every i and fc 

+ (3-38) 

Before proving (3.38) we notice that it implies the required bound for the energy. Indeed, it allows 
us to reduce to the case of one single slope, which was already treated in the first part of the proof. 

It remains to prove (3.38). We first notice that if i,i + k < nt* or i,i + k > nt* , then (3.38) 
follows immediately, since £\ < £2- So we assume that i < nt* < i + k. Then 

"*? = (! + 40 - **) + (1 + li)t* (1 + 

= (£2-£i) l - + {l + £2)^-{£2-li)t* 
>(£2-£i)- + (l + £2)--(£2-£i) f ' ~~ L ' 

n n 

where the last inequality follows from the assumption i < nt* < i + k. The last term is exactly 
(1 + £1)^, and therefore the claim is proved. □ 

4. Uniqueness and convergence of minimizers 
We now prove Theorem 1.2 and Corollary 1.3. 

Proof of Theorem 1.2. To show existence, we take each of the cases k = {1, 2, 3, 4, 5} in turn and 
consider a minimizing sequence [i m for each of them. Note that we only need to prove compactness, 
since each of the T-limits is automatically lower semicontinuous. For E^ 2 \ and E^\ the 

compactness is immediate, since boundedness of E^ k >{ii m ) implies boundedness of f n xd/j, m (x), 
and therefore tightness. For E^ we use the same argument as in the proof of Theorem 3.2 to 
show that J Q x u m (dx) is bounded 

We prove uniqueness by proving strict convexity, either in fj, or in £. Note that convexity of 
£ 1 y E^(^) corresponds to displacement convexity. Also note that the term J Q xfj, = J ^(s)ds 
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is convex in both senses, and therefore we only need to prove strict convexity of the interaction 
terms in each of the E^ k \ 

We treat the cases separately. The functional p i-> E^(p) is strictly convex because the 
function p M> p 2 is strictly convex. Similarly, since s i— > V e s(s) is strictly convex for s > 0, the 
function £ h-» £?( 4 ) (£) is strictly convex. Writing .E^ 2 ) as 

^ (2) (0 = ^ [ 1 V(cW)-Z(s)))dsdt+ f Z{s)ds = c f f V(c(Z(t)-Z(s))) dsdt+ [\(s)ds. 
A Jo Jo Jo Jo Jo Jo 

Since s ^ V(s) is strictly convex for s > 0, it follows that E^ is strictly convex. For E^ a 
similar argument applies. 

The functional E^ is non-strictly convex; but one finds in a straightforward manner that 
p = jC\ , q ^ is the unique minimizing measure. □ 

Proof of Corollary 1.3. Fix the case k G {1, 2, 3, 4, 5}. Given a sequence of minimizers (x™, . . . , x") 
of E^\ we again set p n = ^X)™=i^"- Taking p to be the unique minimizer of E^ given by 
Theorem 1.2, by Theorem 1.1 there exists a recovery sequence /t„ — p along which E^\p n ) 
converges and therefore remains bounded. Since p n are minimizers, E^\p n ) < E^\p n ) also 
remains bounded. By the compactness statement of Theorem 1.1, this imples that p n is compact. 

Therefore p n converges along a subsequence rik to a limit p; the minimality of p nk transfers to 
p, so that p is also a minimizer of the limit E^ (see e.g. [14, Corollary 7.20]). By Theorem 1.2 
this minimizer is unique, p = p, and the whole sequence p n converges. □ 

Appendix A 

A.l. V as a convolution. In this section we prove that V can be written as a convolution. Our 
definition of the Fourier transform and its inverse is 

/(0 := / e-^f(x)d Xl /» := / e 2)ria */(0 <£■ 
Jr Jb. 

We claim that V = W * W, where W := U, U := Vv, and V > 0. 

First of all we notice that, if the function VF is well-defined, then V — W *W. Indeed, since 

v = w * w <s> v = wTw, 

by elementary properties of the Fourier transform and by the definition of W and U we have 

v = w^w = w 2 = u 2 , 

which proves the claim. 

Now we prove that the function W is well-defined. 

We first observe that since V is even, V is real-valued. Moreover, the Fourier transform of V 
can be computed explicitly, as we now show. 

A. 1.1. Fourier transform ofV. We start computing the Fourier transform of the function tp(x) = 
sinh^Trx) = — V( x )- We can rewrite the function ip as 

ip(x) = — — coth(7rx). (A.l) 
7r dx 

In order to compute the Fourier transform of ip, we first note that ip has a non-integrable singularity 
at the origin; therefore <p should be considered a tempered distribution, whose effect on a Schwartz 
function ip is defined as the principle-value integral 

(<p,ip) := -PV [ -(^-coth(nx))ip(x)dx = -PV [ coth(nx)-^(xip(x)) dx. 
J R n\dx ) 7T J R dx 

Below we show that 

coth(7r-)(0 = -icoth(7r£), (A.2) 
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from which it follows by the properties of Fourier transforms that 

m = V(0 = ^^(^ coth «))- 

From the explicit expression of V we can see that it is nonnegative. Since V is an even function, 
it is sufficient to check that it is positive for £ > 0. Writing more explicit the last term in (A. 2) 
we have 

9f (£ coth(7rf )] = coth(7rf) ^ — - = , - . — [ cosh(7rf) . '- f tt | . 

H s v S V v w sinh 2 «) smh«) V smh«)y 

And, since sinhi > t for t > 0, from the previous expression we get 

d 6 U coth«)) > — ^ (cosh«) - 1) , (A.3) 
\ / smn(7r4 ) 

which is positive for £ > since cosht > 1 for every t. Then this proves that V"(£) > for £ > 0. 

This proves in turn that the function U is well defined, and it is even, since V is even. Hence 
the function W is real (once again, since W is defined as the inverse Fourier transform of U and 
U has a non-integrable singularity at the origin, U should be considered a tempered distribution). 

To prove (A. 2) it is convenient to use the representation of the hyperbolic cotangent in terms 
of a series, i.e., 

1 1 ' 



k 2 n 2 + x 2 x ^k 2 TT 2 +x 2 ' 

k— — oo k—1 



then we have 



1 

coth(7rx) = h 2ir ^ 



irx ^ k 2 ir 2 + ir 2 x 2 

k=l 



Therefore, using the previous expression, we have 



coth(7r-)(0= / dx + 27rV / e - 2 ^— — dx. (A.4) 

JM ™ f^J M k 2 TT 2 +TT 2 X 2 



We compute the two integrals in the right-hand side of the previous expression separately. For 
the first term we have 

dx — — isgn(£). (A. 5) 

For the second term we claim that 



/ 



f 



,2 2 T 2 2 dX =~ e ~ 2MC ' S § n (0- (A.6) 
k 7T + -K 2 X Z TT 

To prove this we need some preliminary steps. First of all, elementary calculations show that 

j^h^ * - /, dx - <A - 7) 

Moreover, the Fourier transform of the function / defined as 

tf \ 2a 



a 2 + (2ttx) 2 

is /(£) = e "l^l ; using this formula in our case leads to 

1 . ! 



/ 



e- 2 *** ; 2 2 1 2 2 dx = 4 e - 2 - fc l«l . (A.8) 
k 2 TT z + tt 2 x 2 irk 



Therefore, combining (A. 7) and (A.8) we have 

/ ^ X ,2 2 I 2 2 dx = ^ m = - ^(0 e~ 2 ^, 

J R k 7r + tt z x z 2n nk n 
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which proves the claim (A. 6). Finally, from relations (A. 4), (A. 5) and (A. 6) we have 

oo 

coth]^-) (0 = -i sgn(0 - 2 i sgn(0 ^ e - 27rfc l4l . (A.9) 

fc=i 

At this point we make use of the expression of the hyperbolic cotangent in terms of an infinite 
series, for negative values of its argument, i.e., 

oo 

cothx = -1 - 2^e 2kx , x<0; 

k=l 

then (A.9) reduces simply to 

coth^-XO = isgn(£)coth(7r(-|£|)) = -icoth«). (A.10) 
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